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Many statistical mechanics problems can be framed in terms of random curves; we consider a 
class of three-dimensional loop models that are prototypes for such ensembles. The models show 
transitions between phases with infinite loops and short-loop phases. We map them to CP n_1 
sigma models, where n is the loop fugacity. Using Monte Carlo simulations, we find continuous 
transitions for n = 1,2,3, and first order transitions for n > 5. The results are relevant to line de- 
fects in random media, as well as to Anderson localization and (2+l)-dimensional quantum magnets. 
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Loop models - statistical mechanics problems whose 
degrees of freedom are loops or random walks - are closely 
tied to field theory and more conventional statistical me- 
chanics models, and give an alternative view on critical 
phenomena that has yielded new theoretical approaches 
and tools. More concretely, loops appear as topological 
defects, such as domain walls in two dimensions (2D) and 
vortices in three; as polymers; in the high-temperature 
expansions of lattice models, and in Monte-Carlo ap- 
proaches to quantum problems. Questions about their 
statistics crop up in areas as diverse as Anderson localiza- 
tion [H [2] , turbulence [3] , quantum chaos [4] , cosmology 
[5], optics [5], and frustrated magnetism [7]. 

While there has been great progress in understanding 
2D loop ensembles, the situation in 3D is less clear. This 
is not solely due to the unavailability of exact results: 
many qualitative questions are unanswered. Consider 
as an example random curves appearing in disordered 
media - 'deterministic walks in a random environment' 
(DWRE) . The broad applicability of results for percola- 
tion cluster boundaries to 2D DWRE such as level lines 
of random height functions is well known, as are var- 
ious continuum approaches to this problem [8] 02 [14] . 
But analogous 3D problems, such as the statistics of vor- 
tex lines in random fields, are not as well understood. 
Numerous problems of this kind have been simulated 
El [10], but it has not been clear which are in the 
same universality class, or what the relevant field theo- 
ries should be. 

In this paper we consider a family of three-dimensional 
loop models which are interesting from several points of 
view. A special case has appeared in the study of Ander- 
son localization in 3D [JJ [3] ITTrfTB"] , and results also ap- 
ply to DWRE such as vortices. The loop models are also 
related to (2+l)D quantum magnets. They show transi- 
tions between short-loop phases and Brownian phases in 
which walks can escape to infinity. At critical points the 
loops have a nontrivial fractal structure. 

Our aim is to understand the continuum descriptions 
of these models. We give an analytical mapping to lat- 
tice problems with more conventional (local) degrees of 



freedom. Coarse-graining then yields (compact) CP" -1 
models or supersymmctric variants, field theories which 
in two dimensions have been related to loop models by 
Read and Saleur [14]; see also [T5HT7] . In addition, we 
perform Monte-Carlo simulations of the loop models, ob- 
taining their phase diagrams and accurate values for crit- 
ical exponents. These results support the identification 
of the continuum theory, since we find the expected expo- 
nents in the case (the CP 1 model) where they are known. 
For the CP 2 model we find exponents apparently for the 
first time. 

An important distinction is between oriented and un- 
oriented loop ensembles. While we focus mainly on the 
former, we argue that un-oriented loops are described by 
jjpn-i models. Separately, an extension of the present 
work shows the general applicability of the CP™ -1 and 
RP 11 ^ 1 sigma models in the limit n — > 1 to problems of 
random curves in 3D disordered media |18| . 

Loop models. We consider loop models defined on four- 
coordinated, directed lattices, with two directed links en- 
tering and two leaving each node. A configuration C of 
completely packed, oriented loops is generated by pair- 
ing up the incoming and outgoing links at each node in 
one of the two ways compatible with their orientations 
(Fig. [I]). At each node, one of these pairings occurs with 
Boltzmann weight p and the other weight 1 — p, with 
< p < 1; the assignments are part of the definition of 
the model, along with the choice of lattice. We also give 
the loops a fugacity n. Let \C\ be the total number of 
loops, and N p the number of nodes where the weight -p 
pairing is followed. The partition function is 

Sloops =5> w '(l-P) JVl -'n |C| - (1) 
c 

When the fugacity is a positive integer, it can be re- 
produced by a sum over n 'colours' for each loop. The 
models with n = 1 have the property Z\ oops = 1: In 
this case, the node configurations are independent ran- 
dom variables, and the walks are DWRE. In general, as 
n does not flow under renormalization |19j . we think of 
it as labelling different models, and of p as a parameter. 



This general recipe can be used to construct various 
models, depending on the lattice and node assignments. 
Here we are interested in models with transitions between 
localized and extended phases, and in the universal be- 
haviour at the transition and in the extended phase. Our 
simulations use Cardy's '3D L-lattice' [IT] , which has cu- 
bic symmetry, and a variant, the '3D K-lattice', which 
differs in its link orientations and phase diagram (Fig. [2]). 

Lattice magnet. We rewrite Z\ oops in terms of lo- 
cal 'magnetic' degrees of freedom which can be coarse- 
grained in a fairly straightforward way. In this we are 
inspired by the well-known 0(n) loop models [20] - here 
we obtain instead a lattice CP" -1 model. 

Introduce complex unit vectors z\ = (zf, zf ) on the 
links I of the lattice, and denote the integral over these 
degrees of freedom by Tr (normalized so Trl = 1). Now 
consider a Boltzmann weight which is a product of terms, 
one for each node of the lattice. Labelling the incoming 
and outgoing links at a given node as in Fig. [TJ 

Z = Tr [] (p(4*)(4%') + (1 - P)(4*')(4*)) • 

nodes 

(2) 

This partition function reproduces the sum over loops 
with the right weights. To see this, note that the terms 
in the expansion of the product over nodes are in corre- 
spondence with loop configurations C: 

Z = Tr $>">(1 - pfi-* [] tr H(z lZ }). (3) 

c cec lec 

Here £ is a loop, the outer product zizj is annxn matrix, 
and 'tr' is a trace in this space (the ordering of the last 
product is given by the sequence of links on C). Now, 
since Tr zizj = we are left with one n x n trace, i.e. 
one 'colour' index to sum, per loop. Let Ni be the total 
number of links on the lattice. Then 

Z = ^-X E ^(1-^ = ^^,(4) 

C loop colours 

The Boltzmann weight ^ defines a classical magnet for 
the 'spins' z. In addition to the unitary global symmetry, 
it has the local U(l) symmetry zi — » e^'zi, so our spins 
live not on the sphere \z\ 2 = 1 but on complex projective 
space, CP -11-1 . This space degenerates to a point when 
n = 1, leaving no degrees of freedom. Thus we must 
either resort to a replica-like limit n — > 1, or generalize 
Q to a supersymmetric theory by replacing z with a 




FIG. 1. Pairings at a node (with associated weights), and the 
labelling of links used in Q. 




FIG. 2. Loops on the 3D K-lattice at p = 0. At p = 1, they 
become infinite straight trajectories (crossing at nodes). 

unit supervector of n + k bosonic and k fcrmionic com- 
ponents, tp — (z 1 , z n+k , x 1 , X k )- A nonzero number 
k of fermions leaves the partition function and its loop 
representation unchanged (using Tr-0?/^ = the loop 
expansion goes through as before, with tr — > str) but 
yields more operators, and is necessary to give a nontriv- 
ial theory when n = 1 (or n < 1). 

Field theory. The naive continuum limit of ([2| is the 
(jpn-i moc i c i In a sigma model formulation, with an 
auxiliary gauge field A to remove the unwanted phase 
degree of freedom, the Lagrangian density is: 

£= ^\(d-iA)z\ 2 , with \z\ 2 = l. (5) 

The SUSY version, the CP n+k -^ k model, is got by z -> 
tj). A crucial point in any formulation is that the gauge 
field is compact: the set of gauge transformations z — ¥ 
e^z, A — >• A + d(f> is larger than in noncompact U(l) 
gauge theory as <fi can jump by 27r. This implies that 
Dirac strings of flux 2-7r incur no cost in action, and that 
in integrating over A we must include Dirac monopole 
configurations with quantized charge [21) . 

Work on deconfined criticality [22] has made clear that 
compactness is a subtle issue, so it is useful to have an- 
other route to the continuum limit for the loop models. 
We use the transfer matrix to extract a (2 + l)D quantum 
SU(n) antiferromagnet [23l [25] as an intermediate step 
- for an analogue in 2D, see [13] [TS]. This procedure, to 
be described in [23] > clarifies the compactness of A. 

Finally, an alternative to ^ is to use explicitly gauge- 
invariant degrees of freedom. The two-colour case n = 2 
reduces simply to the 0(3) (classical Heisenberg) model 
via = z^a^z, with a a Pauli matrix, and indeed the 
loop models with n = 2 show the usual 0(3) exponents 
as described below. For general n (without fermions) we 
can use the traceless matrix Q = zz^ — 1/n. 

Correlators. In 3D, the CP n_1 model has a transition 
between a disordered phase and an ordered phase with 
2(n — 1) Goldstone modes (or 2[n + k — 1] bosonic Gold- 
stone modes and k complex fermions). Translating cor- 
relators of gauge invariant operators into loop language 
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shows that the former corresponds to the localized and 
the latter to the extended phase of the loop model. 

Consider Gjv(r), the probability that two small regions 
separated by a distance r are connected by N distinct 
strands of loop. In the localized phase and at a critical 
point all loops are finite, and Gjy(r) is non-zero only for 
even N. G2(r), the probability that two distant points 
lie on the same loop, is proportional to (tr Q(0)Q(r)), 
and higher correlation functions G-2m{t) can be written 
as two-point functions of operators (z 1 z 2 ) M . In the lo- 
calized phase Q is massive and G 2 (r) ~ r^e -1 "^. The 
typical loop size £ diverges on approaching the critical 
point at p — p c , as £ ~ \p — p c \~ u , and at criticality 
G%(r) ~ l/r 1+r ', where 77 is the anomalous dimension 
of the gauge invariant spin Q (or the Heisenberg spin S 
when n — 2). A simple scaling argument [HI relates 
r\ to the fractal dimension df of the critical loops, and to 
the exponent r governing the distribution, P(l) ~ I , 
of the length Z of the loop through a given link: 
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Since we expect 77 to be small, df will be close to 5/2. 
Interestingly, the mean field value 77 = does not give the 
'trivial' fractal dimension of two - this is due to confine- 
ment of z into a gauge invariant composite field. Con- 
finement also has an interpretation as a relation between 
different loop ensembles (essentially worldlines of z ver- 
sus worldlines of S) when n = 2. 

An important basic consequence of the CP™ -1 descrip- 
tion is Brownian behaviour in the extended phase, which 
has been observed but not derived in related problems 
[5H3 HOI LT2]. In this phase, contributions from finite 
strands of infinite walks make Gjv(r) non-zero for both 
even and odd N: Gjsr(r) ~ r~ N . These are Brownian 
exponents (each factor of 1/r is just the probability that 
a random walker visits a given site at distance r from its 
origin) and imply a fractal dimension of two. They follow 
from free field theory for the Goldstone modes, allowing 
for the fact that both the properties of long loops and 
symmetry breaking in the CP n model are sensitive to 
boundary conditions. If walks can end on the boundary, 
'infinite' walks are those that do so: they have typical 
length of order L 2 . With periodic BCs, all walks form 
closed loops: though Brownian, the 'infinite' loops have 
typical length of 0(L 3 ), since a random walker trapped 
in a region of linear size L will on average visit 0(L 3 ) 
sites before re-finding his starting point. The probability 
of a given link lying on an infinite loop is proportional to 
the order parameter, so varies as \p — p c \ close to p c . 

Unoriented loops. Models with unoriented loops are 
also interesting, e.g. in relation to polymers and Z2 vor- 
tices. Similar arguments relate them to i?P n_1 sigma 
models, with real spins. In the models considered above, 
allowing node pairings which do not respect the link ori- 
entation - so loops have no well-defined orientation - 
corresponds to a perturbation SC cx — ti Q T Q, favouring 
real Qs and causing a crossover to i?P n_1 behaviour. 
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FIG. 3. Phase diagrams for the K and L-lattices. Continuous 
transitions are indicated by blue dots and single line and first 
order transitions by red dots and double line. 



Deterministic walks in a random environment. A key 
outcome of this work is the general applicability of the 
(jpk\k moc [ e i t oriented loops in short-range correlated 
random media - most notably various kinds of vortices, 
such as optical vortices [6] , cosmic strings [5] , XY vortices 
in the paramagnetic phase |18j , and 'tricords' in tricolour 
percolation |10) . The striking compatibility of the expo- 
nents in [TO] and in the n = 1 loop model [T^] confirms 
that these problems are in the same universality class. 
Derivations will appear separately |18j . 

Numerical results. We use Monte Carlo simulations 
to study the loop model (1) on the 3D K and L-lattices 
introduced above, for integer n in the range 1 < n < 10. 
We take samples of linear size 32 < L < 100 in units 
of the link length, with periodic BCs. The fugacity n is 
introduced via loop colours, and two types of elementary 
Monte Carlo move are employed: either a change in the 
colour of one loop, or a change in the configuration of a 
node whose links all carry the same colour. Typical run 
lengths involve of order 10 5 Monte Carlo steps of each 
type, per loop or node respectively. For details, see [24] . 

Phase diagrams are shown in Fig. [3| The 3D K-lattice 
is constructed to have only short loops at p = but 
infinite ballistic trajectories at p = 1. Simulations show 
a single transition between a localized and an extended 
phase, which is continuous for n < 3, and first order for 
n > 5 in agreement with a simple large n treatment. 
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FIG. 4. Scaling collapse for n w (p, L) on the 3D K-lattice. 
Main panel: n = 2; inset: n = 3. 
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FIG. 5. First order transitions at n — 5: (main panel) the 
jump in n w on the K-lattice and (inset) in n+ on the 3D 
L-lattice. 



The 3D L lattice is symmetric under p — > 1 — p and 
has only short loops at p = or 1. For n < 3 it has an 
extended phase around p = 1/2, separated by continuous 
transitions from localized phases at large and small p. 
For n > 5 it has only localized phases, and a first order 
transition at p = 1/2. Work is in progress to resolve 
behaviour at n — 4 on both lattices. Previous Monte 
Carlo studies of CP" -1 gHl ES] obtained a first order- 
transition at n = 4. 

We present studies of critical behaviour for transitions 
on the K-lattice at n = 2 and 3. Results for the 3D 
L-lattice are consistent with universality when compared 
with the K-lattice at these values of n and with previous 
work [12] at n = 1. We examine two observables. One is 
the average number n w (p, L) of curves spanning the sam- 
ple in a given direction. The other is the susceptibility, 



which can be expressed [23] in terms of the average num- 
ber n{l) of loops of length I, as x(p> L) = L~ 3 l 2 n(l). 

Empirically, the scaling of the winding number, includ- 
ing finite size corrections, is adequately described by the 
form n w (p,L) = f(x)(l + P m (x)L ViTr ), where: x is the 
scaling variable, x ~ L 1 ^(5p + A5p 2 ); Sp = (p — p c ); v is 
the correlation length exponent; y ltx < 0; P m {x) is a poly- 
nomial of order to; and f(x) is constructed using splines. 
Results are shown in Fig. [4] We fit x(Pi L) m a similar 
manner, using the susceptibility exponent 7. Values for 
y- m (—1.0(3) and —0.6(4), at n — 2 and 3 respectively) 
have large uncertainties, but those of v and 7 are much 
more precise as finite size corrections are small. 

For n = 2 we obtain v = 0.708(5) and 7 = 1.39(1), fit- 
ting to over 300 data points. We believe that the consis- 
tency of these values with previous high-precision studies 
of the 3D classical Heisenberg model [y = 0.7112(5) and 
7 = 1.3960(9) [27]) provides compelling support for our 
identification of the loop model with CP™ -1 . For n = 3 
we find v = 0.50(1) and 7 = 1.01(2). We are not aware 
of a previous determination of exponents for CP 2 . 

In contrast, for n > 5 we find clear evidence of first 
order transitions on both lattices, as displayed in Fig. [5] 
For the K-lattice there is a rapid change in n w (p, L) as p 
passes through p c , developing into a step with increasing 
L. For the 3D L-lattice n w (p, L) — > for large L at all p. 
A transition at p = 1/2 between distinct localized phases 
is signaled by a step in n + = (N p )/(N p + Ni— P ), which 
can be viewed as the internal energy density. 
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